Reducing multi-photon rates in pulsed down-conversion by temporal multiplexing 
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We present a simple technique to reduce the emission rate of higher-order photon events from 
pulsed spontaneous parametric down-conversion. The technique uses extra-cavity control over a 
mode locked ultrafast laser to simultaneously increase repetition rate and reduce the energy of 
each pulse from the pump beam. We apply our scheme to a photonic quantum gate, showing 
improvements in the non-classical interference visibility for 2-photon and 4-photon experiments, 
and in the quantum-gate fidelity and entangled state production in the 2-photon case. 



INTRODUCTION 

True single-photon sources are an essential require- 
ment for scalable applications of quantum information 
processing. Ideally such sources should deterministi- 
cally provide a Fourier limited wave-packet, having one 
and only one photon in a well defined spatio-temporal 
mode, and desirably with a high brightness [I]. Addi- 
tionally, and in particular for the scalability of linear 
optics quantum computing [2J, multiple devices should 
produce identical spectral emission to enable heralded 
and indistinguishable single-photons. Unfortunately, it 
is challenging to engineer single-photon emitters meet- 
ing all of these requirements. Solid state solutions, such 
as quantum dots [31 |3J despite demonstrating both trig- 
gered single-photon emission [5] and high collection effi- 
ciencies [6] , have yet to achieve equivalent spectral prop- 
erties across multiple sources. Systems based on color 
centers in diamond [5] , can suffer from low collection 
efficiencies due to isotropic emission and other candidates 
like atoms [9l [10] and molecules [11] [12] have similar is- 
sues with low output coupling and as of yet have not 
demonstrated high brightness. 

One approach to produce single-photons is to em- 
ploy pair sources, where conditional detection of one 
of a pair of photons is used to herald single-photon 
events. The current state-of-the art for heralded 
single-photon sources is spontaneous parametric down- 
conversion (SPDC) [T3HTrjj . Here a pump laser is used 
to create a pair of photons in a nonlinear birefringent 
material. However, creating single-photons in a pure 
state requires careful engineering of the properties of a 
down-conversion source. First, photon pairs from down- 
conversion are naturally entangled in energy and mo- 
mentum, demanding special techniques to tailor their 
joint correlations and avoid heralding single-photons into 
mixed states [T?H2"0] . Second, since the process is spon- 
taneous, there is a probability of emitting more than a 
single photon into the same spatio-temporal mode [2"T] . 
This effect is intensified when strong pump pulses are 
used to drive the down-conversion. These multi-photon, 
or higher-order, emissions have detrimental effects in ap- 
plications such as quantum cryptography, where they 



can compromise the security of quantum key distribu- 
tion [22) : and in linear optical quantum computing, 
where the noise resulting from these higher-order events 
significantly increases the error rates in quantum cir- 
cuits [21] 123] • 

There are several techniques which can minimise multi- 
photon SPDC emissions. The simplest is to run an exper- 
iment at lower pump power, if one can afford the resulting 
decrease in source brightness. If this is not an option, 
several sources at lower power can be multiplexed. In 
the spatial multiplexing scheme suggested in [25] , several 
down-converters are run in parallel as a photon switch- 
yard. Whenever one source produces a trigger photon, 
the signal photon is switched into one common output 
mode. The efficiency for this scheme has been simulated 
in detail in [26} and demonstrated on a small scale in [27j . 
Alternatively, sources can be multiplexed in time. 

In this letter, we implement a passive temporal mul- 
tiplexing scheme. We show that increasing the repeti- 
tion rate of the pump laser while simultaneously lowering 
the energy of each pump pulse decreases the emission of 
multi-photon events. Our technique allows an improve- 
ment of signal-to-noise ratio for heralded single-photon 
sources from pulsed down-conversion, without compro- 
mising the brightness of single-photon emission. Com- 
pared to previous multiplexing approaches, our scheme 
relies on a simple optical arrangement to perform extra- 
cavity control over the pump power and repetition rate. 
As a quality benchmark, we demonstrate the improve- 
ment in the non-classical interference visibility between 
photons emitted by a single or by independent down- 
conversion sources, as well as the overall quality of a 
photonic entangling quantum gate. 



PUMP REPETITION RATE AND 
HIGHER-ORDER TERMS 

We start by showing the effects of increasing the pulse 
repetition rate of the pump on multi-photon emission in 
a single down-conversion source. The interaction of the 
pump beam with the nonlinear medium during each pulse 
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FIG. 1: Experimental scheme, a) Doubling the pump laser repetition rate. A Coherent MIRA 900 HP mode-locked ThSapphire 
laser outputs approximately 3.8 W of 820 nm pulses with a repetition rate of 76 MHz and a pulse length of approximately 
100 fs. This light is frequency doubled via second harmonic generation in a non-linear bismuth borate (BiBO) crystal giving 
1.53 W centred at 410 nm. An optical delay loop, consisting of two beamsplitters and two high- precision mirrors, splits off half 
of the laser light and feeds it back to the pump mode with a 6.6 ns delay which is equal to half the initial separation between 
two pulses (approximately 2 m). Photon pairs are created via spontaneous parametric downconversion in a type-I phase- 
matched /3-barium-borate crystal (BBO), pumped bidirectionally. The photons are sent through interference filters centred 
at 820 nm with a full-width at half-maximum (FWHM) bandwidth of 2.5 nm before being coupled into single mode optical 
fibers. Photon were counted using standard avalanche photo-diodes. The output of each detector is fed into a commercially 
available counting logic with a coincident time window of 3 ns. The 152 MHz source was pumped with a maximum power 
of 1.53 W, with approximately 50% of this power available in the 76 MHz regime. The two photon coincidence brightness 
were 38.4 counts/s/mW and 40.3 counts/s/mW for the 152 MHz and 76 MHz sources respectively, b) The repetition rate is 
increased by introducing additional 50:50 beam splitters for delay loops, where each delay decreases in length by a half with 
respect to the previous one. This scheme can increases the repetition rate R beyond a factor of two with no further overall loss 
in pump power. 



is described by the Hamiltonian 



H 



h.c. 



(1) 



where a\ and b\ are the photon creation operators into 
signal and idler modes a\ and b\ respectively; and £ 
is the overall efficiency parameter, which represents the 
non-linear interaction strength and carries information 
about spectral properties of the pump laser. Further- 
more, £ is linearly proportional to the electric field am- 
plitude of each pump pulse. The output state of the 
down-conversion process can be written as |29j 
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with A=£ r, where r is the interaction time inside the 
down-conversion medium. From Eq. [2] we see that the 
probability of creating n photon pairs per pulse is given 
by 



P(n) = (l-|A| 2 )|A| 2 " 



(3) 



Thus the joint photodetection rate per second for modes 
a\ and b\ using so-called bucket detectors, i.e. photode- 



tectors without photon-number resolution is 

oo 

C comc = i?^(l-(l-7?n 2 P(n), (4) 

n=l 

where R is the repetition rate of the laser and r\ is the 
product of the detector efficiency and the optical effi- 
ciency to include optical losses and optical coupling [33] . 
As pump power per pulse is increased in an effort to 
increase single-photon brightness multi-photon terms in- 
crease more rapidly leading to lower signal-to-noise ra- 
tios. From Eqs. [3] and [4] the signal-to-noise ratio can be 
approximated by the single pair (n=l) emission over the 
double pair (n=2) emission, 



SNR; 



if 



(1- {l-r,W\\\< 



(5) 



If we now halve the power of each pump pulse, such that 
£ — > C/a/2, while simultaneously doubling the repetition 
rate, R — > 2R, the joint photodetection rate becomes 



2R 



OO 

£ 



(1 - (1 - r]) n ) 2 P(n) 



(6) 



Note that the rate of generating just one pair of photons 
per second is not affected, while events n > 2 are reduced 
by a factor of 2 n_1 . In fact, an equivalent argument can 
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be made for an arbitrary multiple increase in repetition 
rate, m, such that the generic formula for this scheme 
becomes, 

C comc (m) = R}_^ ^7— 15 , (7) 

n— 1 

and the signal to noise ratio becomes 

SNR » ^ = . (8) 

(1-(1-^)2|A| 2 

Perhaps more importantly, the same is true for in- 
dependent sources where two passes through the same 
down-conversion crystal (or equivalently two separate 
crystals) are used to produce independent photon pairs, 
see Fig. [l^) . Since the mathematical argument is equiv- 
alent to above we omit it here and direct the interested 
reader to the appendix. 

APPLICATION IN LINEAR OPTICAL 
QUANTUM COMPUTING 

Since most commercial laser systems do not offer the 
feature of a tunable repetition rate, we developed a sim- 
ple extra-cavity arrangement to do this, see Fig. [T] We 
used this scheme to double the repetition rate of the sec- 
ond harmonic light produced by frequency doubling the 
820 nm line of a 76 MHz TkSapphire laser (Coherent 
Mira 900 HP). Using two 50:50 beam splitters placed 
in series we introduce a delay of 6.6 ns, approximately 
half the time lapse between two original laser pulses, 
in one arm of the pulse doubler circuit while simulta- 
neously lowering the peak pulse power. High precision 
mirror mounts are used to steer light into the down- 
conversion crystal to ensure the same phase matching 
conditions from both arms of the pulse doubling circuit. 
Note that the doubler has to be installed after the second- 
harmonic stage, because of the non-linear dependence 
of the second-harmonic process on pump power. Note 
furthermore, that the delay loop does not require active 
phase stabilisation, and that the temporal delay does not 
have to not be an exact multiple of the master laser, pro- 
vided the pulses are separated by at least the timing jitter 
of the detectors. 

Our scheme reduces the total available pump power 
by a factor of 2 due to the probabilistic recombination of 
the delayed and the original pump mode. This tradeoff 
is acceptable since multi-photon experiments relying on 
SPDC in bulk crystals already have to operate at reduced 
power to minimise noise. Modern SPDC sources based 
on periodically poled crystals only require a few hundred 
microwatt of pump power [3D] and would benefit from 
even higher repetition rates. Figure [IJd) shows how our 
scheme can be extended to repetition rates of 4, 8 etc. 
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FIG. 2: Ratio of 4-photon to 2-photon events for varying pho- 
ton source pump power with pump repetition rates of 76 MHz 
(red triangles) and 152 MHz (black circles). The red solid and 
black dashed lines show the theoretical predictions. Since the 
pump beam at 76 MHz is only one arm of the pulse doubling 
circuit, the maximum power available is equal to 50% of the 
total power at 152 MHz. Errors due to Poissonian counting 
statistics are not visible on this scale. 

without additional optical loss except that introduced by 
the beamsplitters. 

To demonstrate the effectiveness of this scheme we 
measured the rate at which higher-order events from 
SPDC occur as a function of pump power. Figure [2] 
shows the ratio between 4-photon and 2-photon events, 
P(n=2)/P(n=l), as a function of the SPDC pump power 
for pumps at 76 MHz and 152 MHz. These results 
were obtained using spatially multiplexed single-photon 
avalanche diodes in order to count the number of photons 
in each down-conversion mode. The results show that 
in both pump regimes the ratio of 4-photon to 2-photon 
events varies linearly with pump pulse power as predicted 
by Eq. [6] However, there is a clear difference between 
the inclination of the two curves due to a decrease in the 
power available per pump pulse for the pump beam at 
152 MHz compared to that at 76 MHz. The calculated 
slopes of both curves have a ratio of 2.12 ± 0.07, which is 
consistent with the fact that emission of single-photons 
per down-conversion mode is not altered in this scheme. 

It is interesting to test this scheme in the practical ap- 
plication of optical quantum information processing. As 
previously discussed in Refs. [23J [21] higher-order events 
from SPDC are the major factor in degrading the quality 
of entangled states produced by linear optical quantum 
gates: non-single photon inputs have a larger impact over 
the total error-per-gate operation than, for example, er- 
rors due to spatial-mode mismatch. 

We built a photonic controlled-phase (CZ) gate to test 
the effectiveness of the pulse-doubling scheme, see Fig. 3 
c). The circuit implementation is based on Hong-Ou- 
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FIG. 3: Schematic of the experimental setup, a) Input photons at modes ai and b\ produced in a single down-conversion 
crystal; b) two independent photons are heralded at modes a\ and ai by coincident detection of photons b\ and hi- c) A 
controlled-phase gate implemented by Hong-Ou-Mandel interference of the input modes at a partial polarised beam splitter 
(PPBS) using two different photon sources. Input photons are launched from single-mode optical fibres into the quantum gate, 
where one input arm is used to control the temporal delay, At, between the two interfering optical modes. The state preparation 
and tomography is implemented using quarter- (QWP) and half-wave plates (HWP) and polarising beam-splitters (PBS). The 
two input optical modes are superposed at a single partially polarizing beam splitter (PPBS) with nominal reflectivities of 
tih=0 for horizontally, and rjv=2/3 for vertically polarized light respectively. Photons are detected using standard avalanche 
photo-diodes (APD). 



Mandel interference of the two optical input modes at 
a partially polarising beam-splitter (PPBS) — a detailed 
description of the gate can be found at [3T] . Ideally the 
PPBS will have reflectivities of r\H =0 and 7^=2/3 for 
horizontally and vertically polarised light respectively. 
Each input photon encodes a polarization qubit in the 
horizontal and vertical V)) basis. Successful op- 

eration of the gate are post-selected by the detection of 
at least one photon in each output mode, which occurs 
non-deterministically with a probability of 1/9. Condi- 
tioned on post-selection the gate acts to induce a non- 
linear phase shift when both input states are vertically 
polarized i.e. \VV) —> —\VV). Furthermore, the gate is 
entangling and produces the maximally entangled state 
\HD) + \ VA) for an input \DD). 

We used three measurements to assess the effective- 
ness of our scheme in a photonic quantum gate: i) the 
quality non-classical interference, ii) the fidelity of bipar- 
tite quantum states and iii) an entanglement measure 
between two qubits. 

When two indistinguishable photons are superposed 
at a beam splitter they will bunch, that is, preferably 
exit the beamsplitter via the same optical output mode. 
Experimentally, as the relative path difference between 
input photons, At, is reduced (and hence their tempo- 
ral indistinguishability is reduced) the bunching effect 
is seen as a drop in coincident photon detection at the 
two output modes. This is known as Hong-Ou-Mandel 
(HOM) [32] interference and its quality is measured by 
the visibility, V=(C dist - C indlst ) /C dist , where C d ist and 
Cindist are the coincidence counts at the output of the 
beam splitter for distinguishable and indistinguishable 
photon inputs respectively. Non-single photon inputs in- 



crease the likelihood of photons being detected in both 
output modes leading to spurious coincidence events. 
Therefore the quality of HOM interference degrades as 
the ratio of multi-photon events to single-photon events 
from SPDC increases, see Eq. [5] In particular, the visi- 
bility of HOM interference is a direct measure of photon- 
number purity once spatio-temporal mode mismatch is 
accounted for [33] , it therefore serves as a benchmark for 
the quality of any linear optics quantum gate. 

In the case of a controlled-phase gate a visibility of 
80% would be observed for a perfect single-photon source 
and an ideal PPBS with reflectivity 77y=2/3. Fig. [| 
shows how our doubling scheme reduces the detrimen- 
tal effect on the visibility of non-classical interference 
while increasing the pump power. We show this effect 
for two scenarios: i) for input photons produced by a 
single down-conversion crystal shown in Fig. [3] a) and ii) 
photons produced by two down-conversion sources (inde- 
pendent photon inputs) shown in Fig. [3] b) . In the latter 
scenario the signal photon from each pass of the crys- 
tal is heralded with the detection of the corresponding 
idler photon. The experimental data is compared to a 
numeric model created with the Matlab quantum optics 
toolbox by Sze M. Tan [34] and associated linear opti- 
cal quantum computing tools written by T. Jennewein, 
see [33]. This model generates photon number states for 
a source derived from the Hamiltonian in Eq. [T] It prop- 
agates the SPDC state through a series of optical compo- 
nents as described in [26] and detects them with count- 
ing devices that act like bucket detectors, i.e. they click 
for photon events n > 1. The theoretical plots in all 
figures were based on this Matlab model, assuming im- 
perfect non-number resolving detectors with a nominal 
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FIG. 4: Experimental data for Hong-Ou-Mandel visibilities, for varying photon source pump power, in a photonic controlled- 
phase gate. The two photon input state \VV) is interfered at a PPBS by reducing the relative path difference between the 
input photons, a) Interference visibilities for dependent photon inputs and b) visibilities for independent photon inputs for 
varying pump powers. Red triangles show results using a pump laser with a repetition rate of 76 MHz and black circles with a 
repetition rate of 152 MHz. The red solid and black dashed lines show the theoretical predictions and the errors are calculated 
using Poissonian counting statistics. 



efficiency of 60% and a measured PPBS reflectivity of 
?7„=0.682± 0.002. 

Lastly we examine the effects of our pulse doubling 
scheme on the quality of entangled states generated by 
a quantum gate. Although the data presented here was 
taken with the dependent photon source an equivalent 
result would be found for the independent case. We char- 
acterized the entangled state generated by the CZ gate 
using quantum state tomography for dependent photon 
inputs produced by a single down-conversion crystal 
[55] , We prepare the initial input state, \DD) and make 
projective measurements on each output photon with 
the over complete set {\H),\V),\D),\A),\R),\L}}, 
where \D)=(\H) + \V))/y/2, \A)=(\H)-\V))/V2, 
\R)=(\H)+i\V)/V2, and \L)=(\H)-i\V)/y/2 giving 
a total of 36 measurements. The input state \DD) gives 
the maximally entangled output state \HD) + \VA) 
and, as such, is most affected by higher-order photon 
emissions as shown in [531 121] • The measured density 
matrix, p, is reconstructed using a maximum likelihood 
algorithm and compared to the ideal state, pideai- We 
chose the measures of state fidelity, given by, 

F = Tr 2 (y/p^p ideal p^) (9) 

and tangle (concurrence squared) as a test for entangled 
state quality. Figure.]^,) shows the results. We observe a 
stark reduction in the rate of state degradation, whilst in- 
creasing source pump power, as we switch from a 76 MHz 
to 152 MHz pump repetition rate. The results show the 
effect for a dependent downconversion source where the 
goal is to reduce the number of n > 2 events per down- 



conversion mode. 

We also characterised the overall gate performance via 
quantum process tomography as detailed in Ref. [35j . 
The process fidelity is calculated by comparing the re- 
sulting process matrix obtained from experiment, X, to 
that of an ideal process for a CZ gate, Xideal- Figure[5]b) 
shows the effect of increasing laser power on process fi- 
delity, defined equivalently to Eq. [9j 

Finally, we simulated Hong-Ou-Mandel interference in- 
side a controlled-phase gate between photons from two 
independent sources. Figure [6] a) shows the effect on 
non-classical interference visibility as a function of both 
photodetection efficiency and repetition rate of the pump 
laser. Notably, we see that increasing the repetition rate 
by a factor of 10 dramatically increases the interference 
quality. We also show separately the variation of the non- 
classical interference visibility as a function of detection 
efficiency, shown in Fig. |6]b), and pump repetition rate, 
shown in Fig. |6]c). While the interference visibility can 
also be increased by employing highly efficient photodc- 
tectors, this technique is less effective and doing so is 
considered technologically more difficult. 

We note that the maximum possible repetition rate 
used to drive a SPDC source is limited to R = 1/At, 
where At is the coincidence time window which, in turn, 
is dominated by the combined electronic jitter of single 
photodetectors and the coincidence counting logic. Com- 
mercial silicon avalanche photon diodes exhibit a timing 
jitter of typically 400 ps, which can be matched by com- 
mercial counting electronics based on field-programmable 
gate arrays (FPGA). An experiment using these detec- 
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FIG. 5: Photonic CZ gate performance for varying source pump powers and repetition rates, a) Photons from a dependent 
SPDC source are prepared in the initial state \DD). The state quality degrades as source pump power, and hence the relative 
number of higher-order terms, increases. This effect is suppressed by doubling the repetition rate of the source pump laser. 
Data was obtained with a pump laser at 76 MHz (triangles) and 152 MHz (circles). The dashed and solid lines show theoretical 
predictions, b) Similarly the process of the entangling operation performed by the gate degrades with source pump power. Red 
triangles show results using a pump laser with a repetition rate of 76 MHz and black circles with a repetition rate of 152 MHz, 
the red solid and black dashed lines show the theoretical predictions. Errors are calculated using Poissonian counting statistics. 



tors can thus in principle resolve between two down- 
conversion events created by laser pulses at a maximum 
repetition rate of 1 GHz, which can be reached with the 
extra-cavity control detailed in this paper. This is well 
worth considering for SPDC experiments relying on the 
widely used 76 MHz laser we used for our work. However, 
it should be pointed out that femtosecond Thsapphirc 
lasers with 500 MHz and even 1 GHz repetition rate 
are now commercially available }36] and are probably a 
more reasonable approach for a newly designed experi- 
ment which does not require abundant pump power. 



CONCLUSION 

In conclusion, we have demonstrated a simple temporal 
multiplexing scheme that reduces the number of higher- 
order photon events from heralded single-photon sources 
based on SPDC. Our technique improves the signal-to- 
noise ratio as a result of reducing multi-photon events, 
without compromising the brightness and quality of the 
desired single-photon states. We also demonstrated an 
improvement in the performance of a linear photonic 
quantum gate using our source. Our technique could 
be integrated with spatially multiplexed down-conversion 
schemes where multiple crystals and optical switches are 
used to herald single photons into purer Fock-states [27] . 
In such a scheme one could reduce the number of spa- 
tially multiplexed down-conversion crystals to achieve a 
desired signal-to-noise ratio. Future improvements in 



single-photon technologies such as linear optics quantum 
computing and quantum communications will require a 
combination of improvements in sources and detection, 
in particular efficient number resolving photon detection. 
In practice, although number resolving detectors reduce 
the chance of spurious counting statistics, because of op- 
tical loss, they do not remove the need to suppress the 
number of higher-order terms from SPDC. 
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APPENDIX 

The Hamiltonian that describes the creation of pho- 
tons for both forward and backward emission from a 
type-I down- conversion source, see Fig. [T^,), can be writ- 
ten as HE], 

H = i^tial b\ + i&ftal b\ + h.c, (10) 

where the £i and £2 represent the overall efficiencies and 
non-linear interaction strengths for the forward and back- 
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FIG. 6: a) Theoretical simulation of non-classical interference visibility in a controlled-phase gate from two independent photon 
sources. The visibility of interference is shown by the color scale and depends on both the detector efficiency and repetition 
rate of the laser. The simulation assumes an input state of \VV) from independent photon sources pumped with 100 % of the 
available pump power and detected with non-number resolving photodetectors. The free parameter in this plot is the optical 
loss which, fitted to the experimental data, is 40%. b) and c) show cross-sections of the simulated data, shown in a), for varying 
detector efficiency and varying pulse repetition rate respectively. The black marker in these plots shows the experimental data 
point taken from Fig|4]b) 



ward emissions, respectively, they are also linearly pro- 
portional to the electric field amplitude of each pulse; d], 

bp with j = {1,2} are the creation operators of the for- 



ward and backward down-conversion modes. From the 
above equation we obtain the following state, 



\*spdc) = V(l-|Ai| 2 )(l-|A 2 | 2 ) £ A™ 1 K,ni) aiihi J2 A " 2 |n 2 ,n 2 ) Q2 b2 (11) 

711= T12— 



with Ai = £it and Ai = £ 2 r. Therefore, the probability 
of creating n\ and n 2 photons from crystal passes 1 and 
2 per pulse is given by 

P(ni,n 2 ) = (1 - |Ax| 2 )(l - |A 2 | 2 ) | A?" | |A 2 "| . (12) 

For independent sources the presence of photons in 
modes a\ and a 2 are heralded upon a detection event in 
modes b\ and £> 2 respectively. Again, using non-number 
resolving detectors with detection efficiency r), the rate 



per second of jointly heralding photons in modes a\ and 
a 2 is given by 

oo oo 

C comc = RJ2J2 (l-(l-v) ni ) 2 ^-^-v) n2 ) 2 P(ni,n 2 ) 

m — 1 n 2—1 

(13) 

Similarly to the previous argument for dependent pho- 
tons, halving the power per pulse while simultaneously 
doubling the repetition rate gives 



m — 1 ri2 — 1 
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